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Cavity-enhanced radiation-pressure coupling of optical and meclianical degrees of freedom gives rise 
to a range of optomeciianical plienomena, in particular providing a route to the quantum regime 
of mesoscopic mechanical oscillators. A prime challenge in cavity optomechanics has however been 
to realize systems which simultaneously maximize optical finesse and mechanical quality. Here 
we demonstrate for the first time independent control over both mechanical and optical degree of 
freedom within one and the same on-chip resonator. The first direct observation of mechanical 
normal mode coupling in a micromechanical system allows for a quantitative understanding of 
mechanical dissipation. Subsequent optimization of the resonator geometry enables intrinsic material 
loss limited mechanical Q-factors, rivalling the best values reported in the high MHz frequency range, 
while simultaneously preserving the resonators' ultra-high optical finesse. Besides manifesting a 
complete understanding of mechanical dissipation in microresonator based optomechanical systems, 
our results provide an ideal setting for cavity optomechanics. 



Over the past years it has become experimentally pos- 
sible to study the coupling of optical and mechanical 
modes via cavity enhanced radiation pressure, which 
gives rise to a diverse set of long anticipated optome- 
chanical phenomenaP] such as radiation pressure driven 
oscillations |i2j| and, as demonstrated in 2006, dynamic 
backaction cooling [3j HI |5j [H |7j. Moreover, this coupling 
can be exploited to perform highly sensitive measure- 
ments of displacement's, 9J which may enable the obser- 
vation of radiation pressure quantum backaction [10 or 
related phenomena. Major goals in the emerging field 
of cavity-optomechanicsjUl [12], such as ground-state 
cooling [131 [13] necessitate high optical finesse and high 
mechanical quality factors at mechanical oscillation fre- 
quencies exceeding the optical cavity's linewidth. While 
recently impressive progress has been made in creating 
experimental settings in which radiation pressure effects 
can be studied [31 [H [5] [TSl [HI HZ], a prime challenge still 
concerns attaining simultaneously high optical finesse 
and high mechanical Q-factors. Nearly all approaches so 
far have combined traditional cm-sized optical elements 
(mirrors) with micro- or nanoscale mechanical oscillators 
which simultaneously act as mirrors [31 [4] \TE\ [16] (or as 
dispersive element 17J). Other approaches have used the 
intrinsic mechanical modes of the optical elements |18j. 
Yet, in general, the required high reflectivity of the micro- 
element limits its dimension to wavelength size and thus 
sets an upper limit to the mechanical frequencies that 
can be achieved. Also, it is exceedingly difficult to at- 
tain high mechanical quality factors while maintaining 
high optical reflectivity as independent control of opti- 
cal and mechanical degrees of freedom is generally not 
possible. Thus, although remarkably high mechanical Q- 
factors at low frequencies have been obtained [T7], this 
and previous approaches f31 [31 [TS] [THl [TB] have so far 
not succeeded in combining mechanical Q-factors com- 
parable to those achieved in the field of NEMS and 
MEMSpH [201 [211 [221 [23] with the best values of op- 
tical finesse [231 ESI [26]. But exactly the combination 



of both is important for applications such as low loss, 
narrow-band "photonic clocks" \T, ^] and indispensable 
for fundamental studies aiming at approaching and de- 
tecting quantized motion in mesoscopic optomechanical 
systems [IIl[ia[lH]. 

Here we show for the first time independent control over 
both optical and mechanical degree of freedom in one 
and the same microscale optomechanical resonator. The 
demonstration of mechanical normal mode coupling [29] 
within a micromechanical device and the concomitant ge- 
ometry dependence of clamping losses allows a quantita- 
tive understanding of mechanical dissipation. We demon- 
strate monolithic spoke-supported silica resonators com- 
prising a toroidal boundary which allows ultra-high op- 
tical finesse (> 10^) [24 rivalling the best values ob- 
tained in Fabry- Perot cavities [3S]. Independent con- 
trol of their mechanical properties leads to strongly re- 
duced clamping losses allowing for unprecedented me- 
chanical quality factors (e.g. 80,000 at 38 MHz) - ri- 
valling the best published values of strained silicon nitride 
nanoresonators[Tni [20] as well as radial contour-mode 
disk resonators [21] at similarly high frequencies which 
allow entering the resolved sideband regime [30! '3T using 
state-of-the-art optical cavities. It is shown that the mea- 
sured Q-factors are only limited by temperature depen- 
dent intrinsic dissipation [351 [33], which can be reduced 
by low temperature operation. Moreover, the observed 
geometry dependent loss mediated by normal mode cou- 
pling may be of relevance across a wide range of micro- 
and nanomechanical oscillators, for which a detailed un- 
derstanding of dissipation is lacking. 



RESULTS 

Starting point of our analysis are toroidal silica 
microcavites [34] that intrinsically combine ultra-high-Q 
optical whispering gallery modes with around 20 me- 
chanical modes in the 0-100 MHz range fS] which are 
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FIG. 1: Observation of mechanical mode coupling, a Mechanical Q-factors (upper panel) and frequencies (lower panel, 
where solid lines denote results of an FEM simulation) of the radial breathing mode (inset) for varying relative undercut 
u = L/R. The Q-factors were found to be remarkably reproducible for six different samples and strongly geometry dependent 
due to intermode coupling, b Q-factors (upper panel) and frequencies (lower panel) of radial breathing mode and a fiexural 
mode of one toroid reveal an avoided crossing confirming that the dispersion lines of both modes do not cross. The mode 
patterns (radial and fiexural modes) hybridize when approaching the coupling region while the corresponding mode patterns 
switch dispersion lines during the avoided crossing. A coupled harmonic oscillator model (solid lines: coupled Q-factors and 
frequencies; dashed lines: bare Q-factors and frequencies) allows an excellent fit to the data. 



observable in an interferometric readout using the struc- 
ture's optical modes (see methods). Recent work has 
shown that the intrinsic coupling of optical and mechani- 
cal modes of toroidal microcavities via radiation-pressure 
can give rise to the effect of dynamical backaction[5], 
which allows realization of narrow bandwidth photonic 
oscillators |27] or photonic RF down-converters [35^ as 
well as radiation pressure cooling of a mechanical 
oscillator^. The radial breathing mode (cf. Fig. la) in 
particular exhibits strong optomechanical coupling, low 
effective mass (« 10^^^ kg) and high frequency enabling 
the first demonstration of resolved sideband cooling of 
a micromechanical oscillator |30j. To understand their 
mechanical Q-factors which have remained completely 
unexplored so far, we first present a study of mechanical 
dissipation in microtoroids and subsequently show how 
the results from this study can be harnessed to devise 
structures with unprecedentedly low dissipation in the 
frequency range above 30 MHz. 



Observation of mechanical normal mode coupling 

In order to assess the contribution of clamping losses, 
the diameter of the silicon pillar holding the silica res- 



onator (cf. Fig. la) is varied. To this end consecutive 
XeF2 etching cycles, undercutting the sihca structure but 
leaving the silica itself unaffected, are applied. Thus, the 
relative undercut u = L/R (R: radius of the cavity, L: 
length of the free standing membrane, cf. inset of Fig. 
la) can be controlled. The dependence of the measured 
Q-factors of the radial breathing mode (RBM) on the rel- 
ative undercut is depicted in Fig. la. Intuitively, it may 
be expected that higher Q-factors are attained for larger 
undercut due to a reduced clamping area. Interestingly, 
however, the Q-factors of six different microresonators of 
similar size show a strongly non-monotonous dependence 
on the relative undercut which contradicts the simple ex- 
pectation of smaller losses for smaller clamping areas. 
Moreover, the behaviour is remarkably reproducible for 
the different samples. A plot of the measured frequen- 
cies along with an FEM simulation (see methods) of the 
toroids' radially symmetric mechanical modes depicted in 
Fig. la shows excellent agreement of measured and sim- 
ulated frequencies. Furthermore, the FEM simulation in- 
dicates an avoided crossing between different mechanical 
modes as first predicted by Mindlin in 1951 29J. A mea- 
surement with a different sample allows directly observ- 
ing this avoided crossing by employing highly sensitive 
quantum limited optical displacement sensing (see meth- 
ods) which allows to also monitor the fiexural modes ex- 



3 



O 
o 



o 





20,000 






5,000 


a 

■ g15,000 

510,000 
s 

5,000 






4,000 




l-<H 


3,000 
2,000 


* 


0.1 1 10 1D0 1000 ' 
^ ^ pressure (mbar) / 


1,000 











_ 





2,000 angular frequencies fii/f = fij. 



0.75 0.8 0.85 0.9 

Relative Undercut 



0.95 



h « — 




^ 1- i : 



(1) 



Qr/f — Qr/I 



(0) 



r^lyl u and Q-factors 



on the undercut u, these can be 



1,000 



500 



FIG. 2: Linear relation between D and measured Q- 
factors. Mechanical Q-factors (dots with error bars, left 
axis) and corresponding simulated D values (dashed lines, 
right axis) for the radial breathing and a flexural mode of an 
on-chip toroid show a linear relation in the parameter range 
depicted where Q ~ 3 D. Inset: Mechanical Q-factors de- 
pending on the background gas pressure. Below 1 mbar the 
gas damping contribution becomes negligible. The black and 
red lines correspond to fits in the molecular (Q oc p~^) and 
viscous (Q cx p~^/^) gas damping regime. 



hibiting weaker optomechanical coupling. The frequen- 
cies and the corresponding Q-factors of both the radial 
breathing mode and a flexural mode obtained in this sec- 
ond measurement are shown in Fig. [TJd. The frequen- 
cies of both modes decrease for larger relative undercut, 
where the flexural mode shows a steeper slope. Thus, it 
can be observed that the eigenfrequencies of both modes 
approach each other but, as predicted by FEM simula- 
tions, they do not cross. Moreover, both mode patterns 
hybridize and eventually swap their dispersion lines dur- 
ing the avoided crossing (inset of Fig. [T]d). It is exactly 
in the region in which both modes couple, where the Q- 
factor of the RBM is found to be strongly reduced. 
This distinctive behaviour can be well explained with a 
generic model of two coupled harmonic oscillators x^/f, 
representing the bare radial (r) and flexural (f) mode, 
with angular frequencies ^^/i f^nd damping rates T^/[. 
An offset of the silica torus from the equatorial plane of 
the silica disk[27] as well as the single-sided clamping of 
the disk to the silicon pillar gives rise to an apprecia- 
ble coupling rate g between both modes. The coupled 
eigenvalues of the system then read 



1,500 used to asymptotically fit the data (cf. Fig. lb). Using 
the values obtained from this first fit, only g remains as 
Q free parameter in equation ([ij . A least square fit to both 
the measured frequencies and Q-factors yields a coupling 
rate of g ~ 14MHz/27r, much larger than the damping 
rates r,./f/27r. The corresponding coupled frequencies 
fl±/2Tr and Q-factors Q± match the data remarkably 
well (cf. Fig. lb) confirming that indeed both modes 
behave as two coupled harmonic oscillators giving rise 
to an avoided crossing between the RBM and the low-Q 
fiexural mode. Moreover, the data presented here, to 
the authors' best knowledge, manifest the first direct 
observation of normal mode coupling of two different 
mechanical modes within a micromechanical resonator. 



and thus the coupled eigenfrequencies fl±/2TT — 
Re(A±) /2tt and Q-factors Q± = Re(A±) /2Im(A±) can 
be extracted. Assuming a linear dependence of the bare 



Quantifying the clamping losses 

For the Q-factor of the RBM the avoided crossings are 
detrimental. The hybridization with the low-Q flexu- 
ral modes strongly reduces the mechanical quality of the 
RBM in the coupling region due to enhanced displace- 
ment amplitudes in the clamping area. Drawing on the 
quantitative power of FEM simulations, it is indeed pos- 
sible to show that the displacement amplitudes at the 
pillar strongly influence the mechanical Q-factor. In or- 
der to study-and elucidate-the role of clamping losses 
the parameter D, deflned as: 
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is introduced, where i^mech is the total energy stored 
in the oscillator, p denotes the density of silica, and c 
the speed of sound in silica, fim/ZTr is the mechanical 
oscillation frequency and Az (r) denotes the out of plane 
oscillation amplitude at position r, with the integration 
extending over the clamping area Ap. D is the mechani- 
cal quality factor that is expected if the clamping area 
Ap is modelled as a membrane, radiating acoustic energy 

This acous- 
-1 



with a power of P ^ cp il^^ \Az (r) | dA. 



(f2„i!„ech) 



as 



tic loss then leads to a Q-factor of 13 — 1 fp- g 

given by equation ^ . Note that a similar expression can 
be derived following a phonon tunneling approach |36j. 

In Fig. 2 the Q-factors measured for two different 
modes of one microresonator are shown together with 
the simulated values of D. Strikingly, D is directly 
proportional to the measured Q-factors, where Q ^ 3- D, 
revealing good agreement of measured mechanical loss 
and FEM simulations. This unambiguously shows that 
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FIG. 3: Novel optomechanical spoke-supported resonators, a SEM images of on-chip optomechanical resonators 
consisting of a toroidal boundary supported by four spokes (scale bar: 5 fim in each panel). The optomechanical system (outer 
torus) is connected to the silicon chip via silica bridges, decoupling its radial motion from the clamping area in the center of the 
disk. For the depicted design mechanical Q-factors up to 80,000 are obtained at frequencies exceeding 30 MHz. b-h Fabrication 
process of ultra low dissipation optomechanical resonators on a chip. First the spokes pattern is defined in the oxide layer of an 
oxidized silicon chip by UV lithography and HF wet-etching (b). A protective layer of photoresist (c) covers the spokes while 
the structure is under-etched from the sides using an isotropic silicon dry etch (d). The photoresist is removed (e) and the laser 
refiow (f) creates the ultra-smooth toroid microcavity (g) while the spokes remain undamaged as their contact to the silicon 
pillar provides a sufficient heat sink. A final etching frees the spokes and creates the centered silicon support of the structure 
without affecting the ultra- high optical finesse (h). 



the mechanical quality factors are clamping loss limited. 
The prefactor of 3 may be attributed to impedance 
matching conditions between the silica disk and the 
silicon pillar which are not taken into account in this 
model. Note that these, and all following measurements, 
were performed in a vacuum chamber (p < 1 mbar) to 
fully eliminate the possible influence of gas damping (cf. 
Fig. 2). Besides explaining the origin of the mechanical 
dissipation of different modes and proving that plainly 
reducing the clamping area to a minimum is not suffi- 
cient for optimized Q-factors the parameter D moreover 
allows a-priori modelling of the oscillator's mechanical 
quality. This constitutes a powerful, direct tool to design 
high-Q mechanical oscillators with negligible clamping 
losses, which we have used to devise the novel structures 



described below. 



Ultra-low dissipation spoke-supported 
microresonators 



One such geometry is shown in Fig. 3a consisting of 
a toroid resonator with four spokes connecting it to the 
silicon support. Incorporation of the silica spokes into 
a toroidal microresonator is possible by following the 
fabrication steps as shown in Fig. 3b-h. Importantly, the 
design allows independent control of the spokes' width, 
length and position for any major and minor diameter 
of the optical toroid microcavities yielding further 
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FIG. 4: Ultra-low dissipation spokes-resonators, a Measured mechanical Q-factors versus simulated D values for the 
RBM of microresonators with spokes (blue dots) and different modes of high-Q toroids (black dots) at room temperature. 
With the spokes resonators Q-factors of 50,000 can be achieved at frequencies above 20 MHz. The data of the spoke-supported 
samples deviate from the linear behaviour demonstrated in Fig. [5] (indicated by dotted line). Instead they show a room 
temperature saturation (solid line) at Q-factors of Q « 50, 000 which can be attributed to intrinsic losses as discussed in the 
text. The horizontal error bars mark the uncertainty in the calculation of D due to deviations of the actual microresonator 
geometry from the geometry used in the simulation (inferred by SEM imaging). Inset: FEM simulation of the RBM of a 
spokes microresoantor. b Displacement spectrum of a spoke-supported resonator with a mechanical Q-factor of Q = 50, 000 at 
24 MHz. c Transmission spectrum of a spoke-supported microresonator showing optical mode-splitting [37] with an unloaded 
finesse of 230,000. 



degrees of freedom to tune their mechanical frequencies. 
Moreover, using the aforementioned FEM model of 
the clamping loss limited Q-factors, i.e. maximizing 
D, the geometry can be chosen such that the discrete 
eigenfrequencies of the spokes decouple the toroid's 
radial motion from the center of the silica disk, strongly 
mitigating clamping losses. 

Fig. |4]i shows the dependence of the measured Q-factors 
on the simulated values of D for several different samples 
of vastly different geometries with and without spokes. 
Indeed, maximizing D using the spokes support allows 
obtaining reproducibly increased mechanical Q-factors. 
The quality factors which can be achieved at room 
temperature are boosted by approximately a factor of 
ten, e.g. a value of 50,000 is attained at 24 MHz as 
shown in Fig. [4]3. The corresponding clamping loss 
limit as described by our model parameter D can be 
pushed beyond the 10^-level using an optimized design. 
Yet, as shown in Fig. |4^ a saturation of the Q-factors 
becomes evident and the data allow a fit according to 
a saturation model Q^^ = l/{aD) + Q^j. Further 
loss mechanisms which will be addressed below become 
dominant for samples with strongly reduced clamping 
losses and-for room temperature operation-set an upper 
limit of Qsat ~ 50,000. The slope of the linear part 
a obtained from this fit is slightly larger than the one 
obtained from Fig. |2] (dashed line in Fig. |4]i). This 
can be attributed to different impedance matching 
conditions as the sample of Fig. [2] was fabricated using 



a 1 /im silica preform whereas the new samples in Fig. [4] 
consist of a 2 /itm silica layer. 

It is emphasized that the optical properties of the 
toroid microresonators are not affected by the optimized 
mechanical design. The CO2 laser assisted reflow 
process [31] creating the surface tension induced surface 
finish of the silica toroid can still be applied and the 
subsequent XeF2 underetch does not affect the optical 
quality. As a result, the thermal bistability and optical 
mode-splitting [3^, salient features of ultra high-Q 
silica microtoroids |24j . can be observed in the novel 
spoke-supported microresonators. To illustrate this 
point. Fig. [4J: shows the transmission spectrum of a 
spokes resonator with an unloaded optical finesse of 
230,000 displaying mode splitting [37]. At the same 
time, the mechanical design is not adversely affected 
by the laser reflow process, even for very small spokes 
(cf. Fig. 3). Thus, optical and mechanical properties 
of microtoroidal optomechanical resonators can be 
decoupled and optimized independently. The mechanical 
frequencies and Q-factors are controlled by choosing 
appropriate spokes dimensions without affecting the 
choice of minor and major diameter of the toroid defining 
its optical frequencies and finesse [24], 
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FIG. 5: Temperature dependence of mechanical dissi- 
pation. The damping Q^^ of a 38 MHz sample as a func- 
tion of temperature is sliown. Tlie damping decreases from 
= 1/32,000 at room temperature to a record value of 
Q'^ = 1/80, 000 at 410 K. A fit (red line) including losses due 
to TLS (Qxts)' anharmonicity {Q~^iJ and a constant back- 
ground (Q~| ) yields excellent agreement with the data and 
allows extracting an upper bound of Q~j^ — 1/140, 000 for the 
remaining clamping losses. Further suppression of the intrin- 
sic losses could enable Q-factors approaching the extracted 
clamping loss limit of Q = 140,000. Inset: Damping of 
a different sample at 34 MHz illustrating the reliability of the 
fit (red line) of the three aforementioned contributions over a 
wide temperature range from 10 K to 300 K. 

Temperature dependence 

In order to elucidate the apparent room temperature 
limit to the mechanical quality of « 50, 000 and to 
quantify contributions of different origin, we studied the 
temperature dependence of the Q-factors. In Fig. |5]the 
damping of the RBM is depicted as a function of tem- 
perature. The damping of = 1/32,000 measured 
for this sample at 38 MHz and room temperature con- 
tinuously decreases when the temperature is increased. 
In particular, at 410 K a record value of Q = 80,000 is 
attained, on par with the best published values in this 
frequency range [H]. Besides demonstrating a further 
reduction of dissipation, the observed temperature 
dependence allows a detailed understanding of the 
limitations to the mechanical quality. As will be shown 
below, these are due to dissipation caused by two- level 
systems (TLS) common to amorphous materials |32j 
and anharmonicity[33j. The inset in Fig. [s] shows the 
damping Q~^{T) obtained with a 34 MHz sample mea- 
sured from 10 to 300K[3S]. A two parameter fit of the 
TLS losses together with the anharmonic contribution 
(no adjustable parameters) and a constant background 
(clamping losses) yields excellent agreement with the 
data over the full temperature range as well as with 
the parameter range reported in experiments conducted 
with bulk silica|321|3n] (see methods). Using the same fit 
for the high temperature data (295 K to 410 K) displayed 



in Fig. |5] also yields excellent agreement. Moreover, 
the extracted temperature independent background 
sets an upper bound for remaining clamping losses of 
= 1/140,000 experimentally confirming that a 
strong reduction of clamping losses is achieved. For this 
specific sample the D model predicts six times lower 
clamping losses. Thus, it is likely that asymmetries in 
the structure cause additional clamping losses not taken 
into account in the above FEM model leading to the 
observed background of 1/140,000. It is noted that 
thermoelastic damping [40 j for which an upper bound of 
10^^ is found for the resonators considered here, can 
safely be neglected in this analysis. The intrinsic losses 
of silica (i.e. those caused by TLS and anharmonicity) 
are known to exhibit a local minimum at a temperature 
of around 500 K [32] and thus a further increase in 
temperature is expected to yield even higher Q-factors 
than those reported here. 

The quantification of remaining clamping loss limit 
yielding values above 10^ is also very important for 
the prospects of future low temperature experiments, 
as the intrinsic losses in amorphous materials can be 
strongly reduced by cryogenic operation (32j. After 
reaching a maximum at ~ 50 K the intrinsic losses 
decrease again at lower temperatures (cf. data in Fig. 
|4]). It is well known ^2] that the Q-factors reach a 
plateau of Q « 3, 000 before rapidly increasing below a 
drop-off temperature which shows a strong frequency 
dependence. For kHz frequencies, the drop-off temper- 
ature resides in the mK-rangejSj and therefore silica 
is commonly not considered as a suitable material for 
low temperature oscillators. It is noted, however, that 
for high oscillation frequencies the drop-off temperature 
strongly increases and for e.g. 43 MHz (Td = 3K) in 
^32j the intrinsic loss limited Q-factors were measured to 
reach 50,000 at 400 mK while quality factors exceeding 
100,000 may be extrapolated for a temperature of 
300 mK[32j. Thus, the spoke supported resonators are 
very promising optomechanical devices for cryogenic 
operation with the prospect of even improving the 
demonstrated outstanding room temperature Q-factors 
at low temperature, possibly approaching the extracted 
clamping loss limit beyond Q = 10^. 



DISCUSSION 

In summary, we presented novel micron-scale on-chip 
optomechanical resonators that for the first time allow 
independent control of optical and mechanical degrees 
of freedom within one and the same device. Avoided 
crossings of mechanical modes (as first predicted by 
Mindlin[29 ) were observed which enabled a detailed 
understanding of mechanical dissipation, allowing the 
design of optomechanical systems combining mechan- 
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ical Q-factors of 80,000, on par with the best pub- 
lished results in the frequency range above 30 MHz, 
with ultra-high optical finesse. This property yields 
improved performance for applications such as realiz- 
ing narrow band filters and narrow line-width radia- 
tion pressure driven "photonic clocks" [27] . Yet , it is also 
an essential ingredient in studies which seek to demon- 
strate quantum signatures of mesoscopic obiects fTOl l28] . 
To illustrate this point, the ratio of displacement spec- 
tral density provided by radiation pressure quantum 
back-action and the thermal noise spectral density is 
givenlH]: S^yS'^[n,,,] = i+^Z^- A(c/„)XT„,n„, 

= 1 for a linear cavity and ^ = tt for a ring-cavity, 
P: power impinging on the cavity). For cryogenic opera- 
tion at temperature T = 300 mK, a Q-factor of 80,000 
(which seems feasible given our results and the mea- 
surements at 400 mK in [32]) and the realistic param- 
eters A = 1 /xm (wavelength), m = 5ng (effective mass), 
f],„/27r = 20 MHz (mechanical frequency), T = 300,000 
(finesse) and a cavity radius of i? = 35 /im (the cav- 
ity bandwidth K/27r for a toroidal cavity is then given 
by K = c/nRT , where c denotes the speed of light in 
vacuo and n = 1.45 the refractive index of silica) this 
expression reaches unity already at a launched power P 
as low as P = 10 /iW. Thus, radiation pressure quan- 
tum back-action-which so far has remained experimen- 
tally inaccessible-may get within reach. Moreover, at a 
bath temperature of T = 300 mK, and for a mechanical 
Q-factor of 80,000 already one tenth of the radiation- 
pressure cooling rates of Y^j^-K — 1.56 MHz demon- 
strated in silica microtoroids [3D] would enable laser cool- 
ing of the REM in the resolved sideband regime to an 
occupation number of n = ^^nv'^ ^ ^'^^ below unity, 
and would thus allow to venture into the regime of cav- 
ity quantum optomechanics. 

Pertaining to the wider implications of our work, the 
studied aspects of geometry dependent mechanical dis- 
sipation via mode coupling (which is also important for 
nanomechanical beam resonators (42 j) and the demon- 
strated quantification of clamping losses can be of broad 
interest and could directly be applicable to other whisper- 
ing gallery mode resonators made of amorphous [43j .44i 
and crystalline material [251 US] as well as to a wide range 
of other micro- and nanomechanical oscillators, such as 
photonic crystal based systems [46 which are also ex- 
pected to exhibit optomechanical coupling. 



METHODS 
Optomechanical monitoring 

To monitor the mechanical modes' thermal noise, laser 
light is coupled to the high-Q optical whispering gallery 
mode of the silica resonators using a tapered optical fi- 



bre. An external cavity diode laser is locked detuned by 
approximately half the optical line-width k with respect 
to a cavity resonance uiq. Thus, the small, thermal fluctu- 
ations Ax of the cavity radius cause changes in the opti- 
cal path-length which imprint themselves linearly into an 
amplitude modulation AP of the power transmitted by 
the cavity, i.e. AP (x Ax. Spectral analysis of the trans- 
mitted power thus allows extracting the mechanical fre- 
quencies f2iji/27r and dissipation rates T^/2'!t from which 
the corresponding quality factor Q = fl^/T^ is inferred. 
To minimize radiation pressure induced back-action, i.e. 
to rule out any ampliflcation [2] or cooling [S] which would 
lead to a modification of the natural line-width[5j, the 
laser power used in these measurements is kept at a very 
low level, well below the threshold for the parametric 
instability [2]. As an additional check the measurements 
are always performed both on the red and blue detuned 
wing of the optical resonance. The Q-factors from both 
measurements usually agree within less than 5%, indicat- 
ing negligible radiation pressure backaction. 
In order to achieve a higher displacement sensitivity al- 
lowing to observe further mechanical modes of the mi- 
croresonators, and in particular the avoided crossing in 
Fig. an adapted version of the Hansch-Couillaud po- 
larization spectroscopy [4T using a Nd:YAG laser locked 
to cavity resonance is employed |[9j. The polarization of 
the field injected into the taper is chosen such that only 
a small portion couples to a polarization non-degenerate 
whispering gallery mode of the cavity while the larger 
part bypasses the cavity serving as a local oscillator. 
Analyzing the phase-shift of the field exiting the cav- 
ity via a homodyne detection a shot-noise limited sen- 

sitivity [f^m] = ^^^l+4ni/K^ of order 

10^^^m/\/Hz is achieved at the frequencies of interest 
from 10-100 MHz[9j, where P denotes the power imping- 
ing on the cavity. Thus, not only the radial breathing 
mode but more than twenty different mechanical modes 
of the toroid between 0-100 MHz can be monitored. All 
observed modes can be identified with the according 
mode pattern obtained by FEM simulation with a rel- 
ative frequency deviation of less than 5% [S] . 



FEM simulations 

A commercial software package (COMSOL Multi- 
physics) is used for 3D finite element modelling of the 
cavity structures. Using fixed boundary conditions at 
the bottom of the conical silicon pillar, the eigenvalue 
solver provided by the package calculates the eigen- 
frequencies rii„/27r and the corresponding stress and 
strain fields of the meshed structure by diagonalizing its 
coupled linear equations of motion. Summing over all 
elements the total mechanical energy i?moch of the mode 
can be calculated. Typically 100,000 elements are used 
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where in the clamping region the characteristic size of 
an element is of order 200 nm. In order to calculate the 
parameter D, the squared out of plane displacement 
amplitude jA^ (r)|^ is summed over the plane connecting 
the silica resonator and the silicon substrate. 



access to clean-room facilities for micro-fabrication, A. 
Marx for assistance in scanning electron microscopy and 
I. Wilson-Rae for stimulating discussions. 



Intrinsic losses in silica 

For fitting the data in Fig. 5 we employed the known 
model for losses caused by two level systems Qxls 
anharmonic damping The expression for Q^Lg 

originating from microscopic two level systems described 
by a distribution of double- well potentials reads 



1 

TLS 



C-Erf 




dy,(3) 



where Erf(z) = ^ e dx is the error function, T 
denotes the temperature and Q,/2tt the oscillation fre- 
quency. Vq is a scale for the distribution of the potential 
barrier V separating both wells, and Ac scales the asym- 
metry in the depth of the wells (all in temperature units) . 

is a dimensionless exponent, tq a relaxation time and C 
a dimensionless parameter. Measurements spanning from 
11 kHz to above 200 MHz as well as from a few up to a few 
hundred Kelvin allow an excellent fit using the parame- 
ters C = 0.28 ± 0.03, Vb = 667 ± 21 K, Vq/Ac = 7.7 ± 0.7, 
logigTo = -12.2 ±0.18 and C= (1.45 ±0.35) • lO-^iSSj. 
For the fit presented in this work we adopt the above 
mean values for C,, Vq ^nd Vq/Ac while using tq and C 
as fit parameters. The fitted values are logj^pTg — —12.1 
(low temperature fit at 34 MHz), logj^g tq = —12.05 (high 
temperature fit at 38 MHz) and C = 1.8 • 10"^, falHng 
within the parameter range as specified in|33j. The an- 
harmonic contribution Q~^^ is given bv|33|: 



1 ^ ^ C^{T)v{T)T »rth(T) 

T 2pvi{T) i + nv,i(r)' 



(4) 



Here, v-q denotes the Debye and v the sound velocity. 
For the fit in Fig. 5 we adopt uj?, = 0.322 1;^ [32], the 
Griineisen parameter 7^ = 3.6|i33,, the mean lifetime 
of thermal phonons Tth(r) determined in|33j and the 
known values for the specific heat per unit volume 
Cv(r), the sound velocity v{T) and the density p of 
silica (e.g. without any fit parameter. 
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